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Rotation is expected to have an important influence on the structure and the evo- 
lution of stars. However, the mechanisms of angular momentum transport in stars 
remain theoretically uncertain and very complex to take into account in stellar models. 
To achieve a better understanding of these processes, we desperately need observational 
constraints on the internal rotation of stars, which until very recently were restricted 
to the Sun. In this paper, we report the detection of mixed modes — i.e. modes that 
behave both as g modes in the core and as p modes in the envelope — in the spectrum 
of the early red giant KIC7341231, which was observed during one year with the Kepler 
spacecraft. By performing an analysis of the oscillation spectrum of the star, we show 
that its non-radial modes are clearly split by stellar rotation and we are able to deter- 
mine precisely the rotational splittings of 18 modes. We then find a stellar model that 
reproduces very well the observed atmospheric and seismic properties of the star. We 
use this model to perform inversions of the internal rotation profile of the star, which 
enables us to show that the core of the star is rotating at least five times faster than 
the envelope. This will shed new light on the processes of transport of angular momen- 
tum in stars. In particular, this result can be used to place constraints on the angular 
momentum coupling between the core and the envelope of early red giants, which could 
help us discriminate between the theories that have been proposed over the last decades. 

Subject headings: Stars: oscillations - Stars: interiors - Stars: evolution - Stars: indi- 
vidual: KIC7341231 



1. Introduction 



Understanding the effects of rotation on stars is currently one of the key steps to making 
progress in stel lar modeling. Angu lar momentum tra nsport plays a central role in star and planet 
formation (see iBodenheimerl Il995l and iBouvierl 120081 . respectively) . Stellar rotation is intimately 
linked with the star formation process. The distribution of main sequence rotation rates in stars 
indeed arises from the hy drodynamic asse mbly phase and the subsequent interaction of protostars 
and their accretion disks (jShu et al.lll994l ). Rotation also impacts stellar structure and evolution 
by inducing a mixing of the che mical elements inside stars. Meridional circulation i s expected on 



general grounds in rotating stars (lEddingtonlll926l . iMestel 



1953, 



shears can generate hydrodynamic instabilities (see iMaeder fc MevnetlhoOoT 



Mati ns & Zahnll2004r). and inter nal 



Pinsonneault 1997 for 



reviews on mixing in high and low-mass stars, respectively). 

Despite its importance, very little is currently known about the internal rotation of stars, the 
timescales over which it is modified, or the major mechanisms responsible for doing so. The theo- 
retical problem is complex because rotation, convection, and magnetis m are intricately tied to one 



another. Three classes of mecha nisms — hydr odynamic, wave-driven (Charbonnel & Talon 



2005, 



Mathis et allboosl . iMathislbood ^. or magnetic (bough fc Mclntvrel fl998l . ISpruitJll999l . ISpada etal 
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20 ldl ) — could transport angular momentum or induce mixing, and their relative importance is a 
matter of active debate. 

To understand the mechanisms of angular momentum transport in stars better, observational 
data on stellar rotation are needed. Unfortunately, in most cases these are limited to observations of 
surface rotation. Seismology is currently the only tool that allows us to probe the internal rotation 
profiles of stars through the detection of rotational splittings of mode frequencies. Observations of 
solar p-mode splittings have established that the convective envelope of the Sun is in mild differ- 
ential rotation, w hereas the radiative interior rotates as a solid body down to about 0.2 Rq (e.g., 
Schou et al.lll998l ). This discovery revolutionized our understa nding of the solar inte rior and raised 
some questions, such as the thinness of the solar tachochline ([Spiegel Zahnlll992l ). that are still 
open. It also showed that although rotational mixing proces ses have been suc c essful in reproduc 



ing t he observations for massive and intermediate-mass star s (ITalon et al.lll997l. 



2000j) , they predict a solar core rotation rate that is too fast (jPinsonneault et al 



Mae der fc Meynet 



1989 ). As a result , 



other more powerf ul transport processes, suc h as internal gravity waves (|Charbonnel &: Tal on 2005) 
or magnetic fields (IGough Mclntyrdll998l ). have to be invoked. 



However, the Sun alone cannot give us information about the timescale for effective angular 
momentum transport in stellar interiors. Indirect methods, involving the evolution of surface 
rotation in stellar populations, have been the only tools to study this matter so far. The spin down 
of sun-like stars can be used to estimate the timescale over which t he radiative core resp onds to a 
mag netized wind torqu e. Timescale estimates range from ~ 10 Myr (jKeppens et al.lll995l ) to ~ 500 
Myr ( Irwin et al.ll2007 ). with a dependence on mass and rotation rate. The survival of rapid rotation 
in old horizontal branch stars ([Peterson et al.lll983l . lBehrll2003l ) requires that they preserve a core 
reservoir of angular momentum, even when taking into account mass loss on the red giant b ranch. 
Differential rotation in red giants is required to explain this data (pills fc Pinsonneaultll2000i ) . This 
was corroborated by evolutionary models of red-giant stars taking into account self-consistently the 
latest prescriptions of angular momentum transport (jPalacios et al.ll2003l . iPalacios et al.ll2006l ). 



In this paper, we are able to estimate the internal rotation profile of a halo early red gi- 
ant star using seis mology. On e of the main aster oseismic goals of the space missions CoRoT 
(jBaglin et al.ll2006l ) and Kepler ([Borucki et al.ll2010i ) is to probe the rotation profile of stars. These 



miss ions are currently p roviding us with exquisite asteroseismic data (see iMichel &: Baglin 



2012 



and iGilliland et al.ll2010i for overviews of the asteroseismology programs of CoRoT and Kepler, 
respectively). They have alr eady made it poss i ble to extract averaged seismic parame ters of hun- 
dreds of solar-like stars (e.g 



Garcia et al 



2009 



Cha plin et al 



2011 



Verner et al 



2011 



thousands 

of red giants (e.g . iHekker et al.ll2009i . iMosser et al.ll2011bl ). and also of stars in open clusters (e.g. 
Stello et al.l 120101) . With these averaged seismic parameters, the validity of t he scaling relations 
proposed by IKjeldsen fc Bedding! (119951 ) could be observationally established (IHuber et al.1 l201ll ) 
and global stellar parameters such as the mass and radi us were estimated for stars in a wide 
range of evolutio nary stages covering the HR diagram (e.g. IStello et al.ll2009i . iKallinger et al.ll2010i . 



Basu et al.l l201ll ). Moreover, the excellent quality of the photometry provided by these instru- 
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merits enabled us to me asure individual p-mode parameters for many modes (e.g. IChaplin et al 



2010, 


Ballot et al. 


2011), leadin 


(e.g. 


Silva Aeuirre et al. 


2011 





), leadi ng to a very precise modeling of the structure of the observed stars 



rotational splittings once the frequency resolution allows us to do so. 

Evolved stars such as subgiants and red giants are among the most promising objects for 
the purpose of studying core rotation. Indeed, these stars have higher-amplitude modes than 
main-sequence stars because of their higher luminosity, and more importantly, many of their non- 
radial modes are so-called mixed modes. As a star evolves past the end of the main sequence, the 
frequencies of g modes become comparable to those of p modes owing to the high core density. When 
a p mode and a g mode of same degree meet, they are known to avoid each other and exchange 
natures, instead of simply crossing. This phenomenon, known as avoided crossing, is caused by the 
coupling between the p-mode and the g-mode cavities. During this process, both modes have a 
mixed character: they behave as p modes in the envelope and as g modes in the core. These mixed 
modes are very useful because, unlike pure g modes, they have large enough surface amplitudes 
so we can detect them a nd yet are se nsitive to the structure of the core. Mixed modes have been 
theoretically known since lOsakil (119751) discove red them in stellar evolution models . They were first 



obser ved from the ground (Kield sen et al. 



201 Oi ) and Kepler (e.g. ICampante et al 



1995 ) and later from space with CoRoT (jDeheuvels et al 



2011 



probe the structure of the core of subgiants (IDeheuvels Michel 



red giants (jBeck et al 



2011 



evolutionary scenarios (jBedding et al 



Mosser et al 



Mathur et al 



2011 



2011a). These modes enabled us to 



2011 



Metcalfe et al.ll2010l ) and 



201 lal). thus making i t pos sible to discrim inate between 



Mosser et al.l [2012|, see iBeddingi 120121 for a review 



The detection of mixed modes that are split by rotat i on w il l allow us to probe the rotation rate 
of stars ev e n in t heir deepest interior (IKawaler et al.1 Il999i . lLochard et al.1 l2004h . Very recently, 
Beck et al.1 (|2012l ) were able to measure the rotational splittings of mixed modes in three red giant 
stars observed with Kepler and they concluded that the core must rotate at least ten times faster 
than the surface in these objects. 

We here report on the detection of rotationally-split mixed modes in the oscillation spectrum 
of the star KIC7341231, lying on the lower giant branch and observed with the Kepler spacecraft. 
We use these modes to probe the internal rotation profile of this star. The rest of the paper is 
organized as follows: we first give an overview of the atmospheric parameters of the star in Sect. [2j 
We analyze in detail the oscillation spectrum of the star, obtained from one year of Kepler data, in 
Sect. [3l This spectrum is made very complex by the presence of many mixed modes. We show that 
the observed non-radial modes are rotationally split and that the rotational splitting varies from 
one mode to the other, which we interpret as a possible evidence for radial differential rotation in 
the star. To study this hypothesis, we search for a stellar model that is in good agreement with 
both the atmospheric and the seismic constraints of the star in Sect. |U Finally, in Sect. [5j we 
use the observed rotational splittings and our closest stellar models to infer information about the 
rotation of the star. In particular, we perform inversions of the rotational profile using both the 
Regularized Least Squares method and the Optimally Localized Averages technique. 
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2. Atmospheric parameters 



2.1. Existing measurements 

The star KIC 7341231 is also known as HIP92775 or G205-42. It has a V magnitude of 9.96 
(ILaird et al.lll988r). The extremely high proper motion (39.18±0.85 mas/yr in RA and 255.25±1.24 



Latham et alj|2002l ) 



mas/yr in DE. Ivan Leeuwenll2007l ) and radial velocity (—269.16 ± 0.14 kms 1 . 
make it an unambiguous halo star. 

The existing measurements of the star's atmospheric parameters are listed in Table [TJ The 
star was found t o be very metal-po or relative to t he Sun — values quot ed in the literature range 
from —2.18 dex (jLaird et al.1 Il988l ) to —0.79 dex (jAmmons et al.1 120061 ) ■ Because of this, it has 
been included in several catalog studies of low-metallicity stars. However, the atmospheric pa- 
rameters found by these works are often not consistent between each other. For instance, there 
is a wid e spread in the determ ination of th e effective temperatu re. Most values range from 
5470 K (jCasagrande et al.1 l2010l ) to 5483 K (jAmmons et al.l I2006I ). but we also note that two 
studies found a significantly higher effective temperature, around 6000 K ([Latham et al.l 120021 . 
Molenda-Zakowicz et al.ll2008l ). These latter measurements are however dubious because they are 
both found along with a log g around 4.0 for the star, which is completely inconsistent with the 
seismic value of the surface gravity that we derive in the present work (log g = 3.55 ±0.03, see Sect. 
3U). 



2.2. Photometric determination of the star's T c g and [Fe/H] 



Stromgren photometr y is available f o r this star. By using the values of (b — y), mi, and 
c\ that were obtained b y ISchuster et al.1 (120061) for this star along with the Stromgren metal- 
licity scale proposed by iCasagrande et al.1 (120111 ). we derived an estimate of the metallicity of 
[Fe/H] = —1.4 ± .1 dex, confirmi ng that the star is metal-poor. We also combined the Tycho2 



Skrutskie et al. 



BtVt phot ometry (IHog et alJl2000f ) with the infrared JHK$ photometry available from the 2MASS 
catalogue 
scribed by 



20061) an d applied the InfraRed Flux Method (IRFM) method as pre- 



Casagrande et al.l (|2010r ) to derive an estimate of the effective temperature of the star. 



We thus obtained T e g- = 5470 ± 150 K. To determine the error bar on this measurement, we took 
into account the uncertainty on the star's metallicity and surface gravity, although the IRFM is 
only weakly dependent on these quantities. The largest sources of uncertainty come from the red- 
dening (depending on the reddening calibration that we use, we obtain < E(b — y) < 0.02) and 
the errors in the photometric observations. 
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Table 1: Atmospheric parameters of KIC7341231 found in the literature or derived in this study. 



T cff (K) 5470 ± 150 Present study (IRFM, Sect. |22J 

5233 ± 50 [100] tt Present study (HIRES spectroscopy, Sect. [2J3 

5915 ± 182 Molenda-Zakowicz et al. (2008) 

5483 Ammons et al. (2006) 

6000 Latham et al. (2002) 

5300 Cavrel de Strobel et al. (2001) 
5251 Thorburn (1994) 

5301 Laird et al. (1988) 



[Fe/H] (dex) —1.4 ±0.1 Present study (Stromgren calibration, Sect. 12.2 

-1.64 ± 0.05 [0.1]" Present study (HIRES spectroscopy, Sect. E3D 

-1.68 ±0.19 Molenda-Zakowicz et al. ( 2008) 

-0.79 Ammons et al. ('2006 ') 

-1.51 Pilachowski et al. (1993) 

-2.18 Laird et al. (1988) 



log g 3.55 ± 0.03 Present study (seismology, Sect. 13.1 

4.06 ± 0.29 Molenda-Zakowicz et al. (2008) 

3.7 Charbonnel fc Talon (2005) 

4.0 Latham et al. (2002) 



vsini (kms -1 ) < 1.0 ± 0.5 [1]" Present study (HIRES spectroscopy, Sect. I2T31) 
" The errors in brackets are obtained by multiplying the internal errors by a factor two (see text). These 

values were used in the present study. 
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Fig. 1. — Spectrum of KIC7341231 obtained with the HIRES spectrograph at the Keck Observatory 
around the Mg B triplet lines. The red dashed line corresponds to the best-fit model. 
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2.3. A new spectroscopic determination of the atmospheric parameters of the star 



Because of the wide variations between the atmospheric parameters f ound for KIC 734 1 231 by 
previous studies, we re-observed the star using the HIRES spectrograph (jVogt et al.lll994l ) at the 
Keck Observatory on 26 July 2011 in order to perform a dedicated spectroscopic analysis. The 
spectrum we obtained has a SNR of about 100 and the B5 decker was used, yielding a spectral 
resolution, R = 52,000. The analysis of the spec trum was carried out using Spectroscopy Made Eas y 
(SME), a spectral synthesis modeling program (jValenti &; Piskunovlll996l . IValenti &; Fischerll2005l ). 
The 1-D LTE radiative transfer code built into SME was used to fit the continuum-normalized 
spectrum of the star. 

It is well known that the determination of log g from spectroscopy alone is not reliable (e.g. 



Smallevl 120051 ) . Instead, global seismic properties of the star (large separation of p modes and 
frequency of the maximum amplitude) provid e a significantly more prec ise estimate of the surface 
gravity through well-established scali ng laws (Kie dse n fc Bedding 1 1995T ) which have recently been 



further validated from observations (H uber et al 



20 111 ). For our analysis, we fixed the surface 



gravity to the seismic value of log g that we obtain from the Kepler oscillation spectrum in the 
present work (log g = 3.55, see Sect. 13. 1|) . This has the advantage of further constraining the other 
spectroscopic parameters such as the effective temperature or the metallicity. 

We searched for optimal values of T e ff , vs'mi, and [Fe/H], using a Levenberg-Marquardt algo- 
rithm to minimize the xt parameter between the model and observed spectrum. After obtaining 
one set of model parameters, we varied the effective temperature by ±400 K as input to new SME 
trials to explore the parameter degeneracy. In each case, the best-fit model returned to the val- 
ues listed in Table [J_J within half the value of the stated uncertainties. A wavelength segment of 
the spectrum of the star, spanning the Mg B triplet lines, is shown in Fig. [TJ The model fit is 
overplotted as a red dashed line in this plot. 

We note that this star is very metal poor and it is quite possible that systematic sources 
of error exceed our estimates of the formal errors. For example, the stellar atmosphere model 
used to generate the synthetic spectrum was obtained by interpolation from a grid of Kurucz 
model atmospheres, which may not be robust for such a metal poor star. Furthermore, the atomic 
line data (oscillator strengths and broadening coefficients) were originally tuned to match the high 
resolution solar spectrum. To account for these likely sources of systematic errors in our analysis the 
formal internal errors have been multiplied by an arbitrary factor of two to provide the uncertainties 
listed in Table [TJ We thus obtained an effective temperature of T e fj = 5233 ± 100 K, lower than 
but in marginal agreement with our result from photometry. In the following, we chose to use 
the photometric estimate of T e g as a reference instead of the spectroscopic measurement because 
the former is directly tied to the fundamental definition of the effective temperature. However, we 
have checked that all our results remain unchanged when considering the spectroscopic temperature 
instead. 
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3. Seismic properties of KIC7341231 



3.1. Kepler observations 

The star KIC7341231 was observed with the Kepler spacecraft over a period of one year 
(quarter^] Q5, Q6, Q7, and Q8 of Kepler observations) with the short-cadence mode (58.84876 
s), allowing the detect ion of solar-like osc illations. The lightcurve was processed using the Kepler 
mission data pipeline (I Jenkins et al.ll201Cll ) and corrected for outliers, occasional jumps and drifts 
following iGarcia et al.l (|201ll ). We also corrected the lightcurve for long-period instrumental drifts 
by applying a high-pass filter (subtraction of a triangularly smoothed version of the lightcurve over 
a width of 1 day). 

The oscillation spectrum of the s tar (shown in Fig. [2]) was obtained from the time series by 
using the Lomb-Scargle periodogram (|Scarglelll982l ), It exhibits a clear excess of power between 
300 and 600 /xHz, which corresponds to the signature of high-order pressure modes. The frequency 
of the maximum signal is i- Wx = 406 ± 3 uHz ( this value is estimated in Sect. I3.3|) . Using the 



scaling relations suggested bv lBrown et al.l (|1991l ). we can derive an estimate of the surface gravity 
from f max and T c g . Using our photometric estimate of the effective temperature (T e g = 5470 ± 150 
K, see Sect. I2.2H . we obtained a seismic estimate of log g = 3.55 ± 0.03. We note that even when 
considering instead our spectroscopic estimate of the effective temperature (T e g = 5233 ± WOK, 
see Sect. I2.3p . our estimate of logg is not significantly modified (we then obtain log g = 3.54). 

An autocorrelation of the power spectrum provides a first estimate of the mean large separation 
(Ay) ~ 29 /zHz. Fig. [3] shows an echelle diagram of the observed power spectrum folded with this 
value of the large separation. Two clear neighboring ridges stand out (at an abscissa around 7 /uHz 
and 3 /xHz in Fig. [3]), which are readily identified as the I = and I = 2 ridges, respectively. No 
actual I = 1 ridge can be identified in Fig. [3l Instead of lying along a ridge, the I = 1 modes 
appear scattered in the echelle diagram. This is caused by the fact that these modes all have 
a mixe d behavior and therefore no longer follow the asymptotic approximation of high-order p 
modes. iDeheuvels &: Michell (|2010bl ) have shown that I = 1 avoided crossings induce distortions in 
the whole I = 1 ridge, owing to the strong coupling between the p-mode and the g-mode cavities for 
these modes. The existence of numerous I = 1 g modes in the frequency range of the observation 
can therefore account for such a scattered I = 1 ridge. This hypothesis is corroborated by the low 
value of log<7 found for this star, which indicates that it is probably in the post main sequence 
stage. 

We note that in this case, we also expect the I = 2 modes to have a mixed behavior. However, 
I = 2 p modes are known to propagate less deeply inside the star than I = 1 p modes. As a result, 



lr The Kepler spacecraft operates over 4 seasons (quarters) each year. At the beginning of each quarter, data from 
the previous quarter are downlinked to the ground, new target tables are uploaded, and the spacecraft is rotated by 
90 degrees to reposition the solar arrays. 
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Fig. 2. — Power density spectrum of KIC7341231 computed with one year of Kepler observations. 
Left: Spectrum smoothed over a 2-/xHz boxcar. The red curve corresponds to a fit of a function of 
the type B{y) (see text) to the spectrum. The contribution from the photon noise is overplotted 
in gray. Right: Power density spectrum over the frequency range where oscillation modes are 
observed. 




Fig. 3. — Fchelle diagram of KIC7341231 folded with a large separation of Av = 29 /jMz. For more 
clarity, the power spectrum was binned over a 0.2-^/Hz boxcar and clipped at a maximum of 50 
times the noise level. The whit e diamonds correspond to the mode frequencies predicted by the 
method of lBenomar et ajj (|2012l ) (see text). 
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the evanescent zone that couples the p-mode and the g-mode cavities is wider for I = 2 modes and 
the coupling is therefore much weaker. This causes the I = 2 modes to be more efficiently trapped. 
In this case, the modes that are trapped in the core have a much larger inertia, an d thus a much 



longer lifetime, than the modes tra pped in the envelop e (jDziembowski et al.ll200ll ). For stars at 



the bottom of the red giant branch, iDupret et al.l ( 20091 ) found a lifetime of the order of 10 3 to 10 
days for I = 2 modes trapped in the core. Even with one year of Kepler da ta, these modes are thus 
expected to be unresolved. This is important because IDupret et al.1 (J2009J) showed that the heights 
of unresolved modes (i.e. their maximum power spectral densities) are inversely proportional to 
their inertia. As a result, the I = 2 modes that are trapped in the core should have very small 
heights. This explains why we detect only the I = 2 modes that are trapped mainly in the envelope, 
i.e. those that are located near the I = 2 ridge of p modes in the echelle diagram. 

One particularly striking feature of the oscillation spectrum is that the I = 1 modes seem to be 
split in two peaks. This phenomenon can even be seen with the naked eye in Fig. [3j This is clearly 
the signature of stellar rotation with a high inclination angle between the rotation axis and the line 
of sight. It opens very exciting possibilities because many of the non-radial modes have a mixed 
nature. The modes that behave mainly as g modes are more sensitive to the rotation in the core, 
whereas the modes trapped in the acoustic cavity are more sensitive to the rotation in the envelope. 
With KIC7341231, we thus have the opportunity to probe the internal rotation profile and look 
for possible signatures of differential rotation in radius. We note that a detailed observation of 
the lightcurve did not show any clear stellar modulation produced by the crossing of spots at the 
surface of the star, although the inclination angle is expectedly high. This is compatible with an 
old star without any strong magnetic field in the outer layers that could trigger surface activity. 



3.2. Identification of the modes 



First, we need to identify the I = 1 mixed modes in the observed power spectrum. iBenomar et al 



(|2012l ) prescribed a method to i dentify the modes in oscillat ion spectra with avoided crossings by 
using a toy model proposed by iDeheuvels Michell (|2010bl ) . If the p and g-mode cavities were 
uncoupled, their eigen modes would coexist independently and no avoided crossings would occur 
(jAizenman et alJll977l ). We will further denote as n modes the theoretical uncoupled p modes and 
as 7 modes the uncoupled g modes. To interpret the observed spectrum, which can be very complex 
because of mixed modes, it is very convenient to consider the uncoupled modes as harmonic oscilla- 
tors and to introduce a coup ling term between them to simulate the effects of the evanescent zone 
( Deheuvels fc Michell l2010bl ). Indeed, this toy model translates into a simple eigenvalue problem 



that can readily be solved to estimate the eigenfrequencies of the coupled system. These frequencies 
can then be fitted to the observed ones (central points of the observed rotational multiplets), by 
adjusting the frequencies of th e theoretical tt modes, those of the 7 modes and the strength of the 
coupling (jBenomar et alj|2012f ). The convergence of such a method is made easier by tak ing into 
account the fact that the uncoupled modes roughly follow the asymptotic approximation of lTassoul 
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(1980) (it modes are approximately equally spaced in frequency and 7 modes in period). 



In the present case, we found a satisfactory fit to the observations by considering 15 / = 1 7 
modes between 310 and 580 /uHz, with a mean period spacing of (Alii) = 107.1 ±2.3 s, and a set of 
7r modes with a mean large separation of (Az/) = 28.9 ±0.2 /^Hz. The frequencies of the I = 1 modes 
that we obtained using this simple model are overplotted in Fig. [3l We thus confirm that most of 
the modes detected outside of the I = and I = 2 ridges can be identified as I = 1 mixed modes and 
we obtain a list of guessed frequencies for the observed I = 1 modes. In a complementary way, we 
analyzed the mixed m o des u sing the asymtotic relation that was applied to hundreds of red-giant 
stars by iMosser et ajj ([2012]) . We obtained estimate s of the I = 1 mixed mode frequencies that 
are very close to the ones found with the method of iBenomar et al.l (|2012l ). We derived a period 
spacing of 112.8 ± .3 s for the I = 1 mo des, in marginal agreement with the value obtained by 
using the method of lBenomar et al.1 (12012 ). 



Unfortunately, these two methods can not be applied to I = 2 modes. Indeed, the coupling 
between the cavities is too weak and as mentioned above, only the modes that are close to acoustic 
modes are detected. However, knowing tha t the period s pacin g of I = 2 modes is such that 
AII2 = AIIi/\/3, the asymptotic relation of IMosser et al.l ( 20121 ) allow us to predict that I = 2 
avoided crossings should occur around the tt modes of radial orders n=ll, 12, 13 (frequencies 
around 327, 355, 384 /xHz). We come back on this matter in Sect. HI using best-fit models of the 
star. 



3.3. Estimating the rotational splittings of the observed non-radial modes 

Rotation is known to lift the degeneracy between the non-radial modes of same radial order 
n and degree I but different azimuthal order m. In the case where the rotation of the star is 
slow enough so that the effects of the centrifugal force can be neglected, a first-order perturbation 
approximates well the effects of rotation on the mode frequencies. The star we study here is in this 
case because it has a low projected surface velocity (v s'mi < 1 ± 1 kms -1 ) and a high inclination 
angle (this latter point will be confirmed in Sect. I3.3.2|). which places it significantly below the 



limit of validity of the perturbative approach (see e.g. ISuarez et al.l I201Q ) . If we further assume 



that the rotation profile is spherically symmetric, the frequency of the (n, I, m) mode is given by 

Vn,l,m = v n ,l,o + m5u U: i (1) 

where b~v n .i is known as the rotational splitting and can be expressed as a weighted measure of the 
star's rotation rate J7(r) 

r R 

<5^ = (2vr)- 1 / K n>l (r)n(r)dr (2) 
J 

The functions K n ^(r) are the rotational kernels of the modes; they depend on the equilibrium struc- 
ture of the star and on the mode eigenfunctions. The expression of rotational splittings for spher- 



ically symmetric rotating stars, which is given by Eq. [H was first obtained bv lCowling &: Newing 
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(119491) and 



Goupi] 



2011 



Ledquxl ()1949l ) (for a review on the effect of rotation on the mode frequencies, see 



). We note that Eq. Q] implies that the components of a rotational multiplet are 



e.g. 

expected to be uniformly spaced by the splitting 5v n j. 

To estimate the rotational splittings of the modes in the oscillation spectrum of the star, we 
followed and adapted a proced ure that was designed to analyze the oscillation spectra of solar-like 



pulsators observed from space ([Fletcher et al 



2006 



Appourchaux et al.l 120061). which was success 



fully applied to CoRoT targets (e .g. Appourchaux et al. 20081 . Deheuvels et al. 201C) and Kepler 
targets (e.g. ICampante et al.ll201ll ). 



3.3.1. Modeling the Power Spectral Density (PSD) 

The observed power spectrum is distributed around a m ean profile P(y), following the statistics 
of a x 2 with 2 degrees of freedom (jDuvall &: Harvevlll986l ). The profile P(v) can be split into two 
components: 

• the background B(v), which is composed of the photon noise and the contribution from 
granulation. The modeling of B(i>) is described in Appendix lA) 

• the contribution from the stellar pulsations to the PSD, which is described below. 

Solar-like oscillations are stochastically excited by the turbulent motions in the outer connec- 
tive envelope. They are intrinsically damped and their profiles in the PSD can be modeled as a 



Loren tzian function with a linewidth inversely proportional to the mode lifetimes (jDuvall &: Harvey 



19861 ). We can thus write 



pt \ = ST ai,m(i)Hn,i ( o\ 

1 + A[v _ Vnl + m5Vnl] 2/ V 2J W 



n.l.m 



where v n ,h H n ,l an d r Hj ; correspond to the frequency, height, and linewidth of the m = component 
of the (n, I) multipl et. Inside a rotational m ultiplet, the height ratios only depend on the inclination 



angle i of the star (jGizon &; Solankil 120031 ) and correspond to the terms a/ jm (i) in Eq. [3l Unlike 
previous analyses of that type, here we allow the rotational splitting 8v n i to vary from one mode 
to the other. The peaks are indeed narrow enough so we can measure individual splittings, and 
hence obtain observational constraints on the variations in the rotation rate with radius. 

Several assumptions that are commonly made while analyzing the spectra of main-sequence 
solar-like pulsators are no longer expected to be valid in the case of a star as evolved as this one. 
For instance, we usually assume that in each overtone, all the modes share a common linewidth T. 
Here, numerous non-radial modes have a strongly mixed behavior, which affects their inertia and 
therefore their lifetimes. We thus considered the linewidths of all the modes as free parameters of 
our fit. Similarly, for main-sequence stars, the ratio between the height of a non-radial mode and 
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the height of the closest radial mode i s usually fixed to a th eoretical value, determined from the 
stellar limb-darkening profile (see e.g. iGizon &: Solankill2003l ). In the present case, we know that 
some modes are mainly trapped in the core and have longer lifetimes. We can therefore not expect 
all the modes to be resolved, which forbids us to use these theoretical visibility ratios. All the mode 
heights are thus left free in our fit. 



3. 3. 2. Extraction of the mode parameters 
The mode parameters have then been independently determined by seven te ams, either by using 



the maximum likelihood estimation (MLE) method (see lAnderson et al.lll990l for more details on 



how this method is applied to oscillation s pectra) or by perform ing a maximum a posteriori (MAP) 



estimation based on Bayesian priors (e.g. iGaulme et al.ll2009i ). Generally, a simultaneous fit of all 



the modes is preferred because numerous parameters can be considered to be common to several 
modes (e.g. linewidths, height ratios...). However, in our case, most of these assumptions had to 
be abandoned and the only parameter which is shared by the modes is the inclination angle of the 
star. For this reason, several teams chose to perform local fits of the modes to gain in flexibility, 
leaving the inclination angle free for each mode in the fit. 

For / = 1 modes, the fitting teams obtained estimates of all the mode parameters (including 
the rotational splittings) in close agreement with one another. In order to obtain a robust list of 
rotational splittings, we selected among the fitted I = 1 modes a subset of 15 modes for which 
at least six of the seven t eams agreed on all parameters within 1-a error bars, as prescribed by 



Appourchaux et al.1 (120081 ) . For these modes, the fits all converged toward a high inclination angle 
(between 70° and 90°). This was expected because the overall profile of the I = 1 modes (which 
appear to be split in two prominent peaks) can only be explained by i ~ 90°. 

The case of I = 2 modes is more complicated. For a rotation axis perpendicular to the line 
of sight, the m = ±1 components of the I = 2 multiplets vanish and we thus expect to detect 
rotationally-split I = 2 modes as triplets composed of the m = and m = ±2 components. In 
practice, the profiles of / = 2 modes are clearly split by rotation only in the three radial overtones 
around the maximum of the signal. Outside of this region, the signal-to-noise ratio is too low 
to detect the signature of rotation. For these latter modes, there was a poor agreement between 
the different teams concerning the rotational splittings and we thus chose to discard them in the 
following. In two of the three overtones around f maX) the I = 2 modes (around 409 and 439 /iHz) 
exhibit the expected behavior and their rotational splittings were successfully estimated (see Table 

E). 

The third of these overtones (in the neighborhood of the radial mode at a frequency of 384.5 
/iHz) shows some peculiar features. The I = 2 mode in this overtone is expected to be located 
between the radial mode and an I = 1 mixed mode at a frequency of 379 /iHz. Fig. 0] shows that 
at least five peaks can be detected in this region instead of three, as would be expected. Besides, 
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when fitting these peaks as the components of one single I = 2 mode, we find a rotational splitting 
of about 0.6 fiHz, which is much too large compared to the splittings of the two other I = 2 modes 
(0.13 and 0.18 //Hz). The most reasonable explanation is that we are in fact detecting an I = 2 
avoided crossing and that two I = 2 modes are present (one around 380.5 /iHz and another one 
around 382.5 /iHz) instead of one. This hypothesis is plausible since the use of asymptotic relations 
had predicted an I = 2 avoided crossing around this frequency (see Sect. I3.2f) and it will be con- 
firmed by the modeling of the star in Sect. HI But even in this case, the multiplets do not have the 
expected pattern. The frequencies of the components in the multiplets appear to be asymmetrical, 
contrary to what would be expected for such a slow rotation. There seems to be an asymmetry also 
in the amplitudes of the components. For these reasons, we failed to obtain rotational splittings for 
these two I = 2 modes in avoided crossing. A theoretical study of the interaction between rotation 
and avoided crossings for I = 2 modes is under way and could possibly explain the observed profile 
of the multiplets. However, at this point, we preferred to discard these two modes. 

A last global fit was performed, using only the modes for which at least six of the seven fitting 
teams agreed to estimate the inclination angle of the star. The idea was to remove the possible 
influence of spurious non-radial modes when determining the angle. We thus obtained an inclination 
angle of i = 85 ± 5° and a robust, reliable list of rotational splittings (given in Table [2]) that can be 
safely used to derive information about the star's rotation profile. The obtained splittings range 
from 0.13 to 0.41 /xHz, with error bars of 0.03 //Hz on average. It is therefore clear that there are 
significant variations of the rotational splitting from one mode to another. To illustrate this, Fig. 
[5] shows the profiles of four I = 1 multiplets that have different rotational splittings. The profiles 
of the two I = 2 multiplets for which the rotational splitting could be estimated are plotted in 
Fig. EJ This shows that the interior of the star rotates differentially in radius. To investigate this 
further, it is necessary to find a stellar model in order to establish the relation between the obtained 
rotational splittings and the trapping of the modes. 



4. Seismic modeling of KIC7341231 

To obtain information about f2(r), we need to have access to the rotational kernels of the 
modes and therefore to a model of the star. Our aim here is to find a model offering a reasonable 
agreement with both the atmospheric constraints of the star and the observed mode frequencies, 
so that we can use the rotational kernels of this model to estimate the rotation profile of the star. 



4.1. Modeling red giant stars using the frequencies of mixed modes 



It is known that the modeling of stars with avoided crossings is co mplex because the timescal e 
of these phenomena is very short compared to the evolution timescale (jDeheuvels &: Michelll2010al ). 



-15- 



Table 2: Estimates of the frequencies and rotational splittings of the detected modes for 
KIC7341231, obtained by fitting Lorentzian functions to the mode profiles. The rotational split- 
tings are given only for the modes for which at least six of the seven teams agreed to within 1-cr. 



1 


v n j (//Hz) 


5v n j (//Hz) 





271.150 ±0.072 


n.a.* 





299.057 ±0.101 


n.a. 





327.239 ±0.051 


n.a. 





355.869 ± 0.040 


n.a. 





384.498 ± 0.030 


n.a. 





413.478 ± 0.035 


n.a. 





442.596 ± 0.040 


n.a. 





472.016 ±0.065 


n.a. 





501.328 ±0.068 


n.a. 





531.517 ±0.202 


n.a. 
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1 


392.160 ± 0.019 


0.252 ± 0.019 


1 


400.559 ± 0.022 


0.156 ±0.024 


1 


412.584 ±0.012 


0.323 ±0.013 


1 


425.766 ± 0.038 


0.216 ± 0.040 


1 


434.839 ± 0.023 


0.262 ± 0.022 


1 


450.966 ± 0.034 


0.292 ± 0.032 


1 


460.377 ± 0.036 


0.270 ± 0.036 


1 


476.696 ± 0.072 


0.285 ± 0.063 


1 


488.244 ± 0.044 


0.243 ± 0.045 


1 


517.382 ± 0.058 


- 


1 


534 012 ± 096 




1 


547.940 ±0.187 




2 


295.212 ±0.300 




2 


323.649 ±0.218 




2 


352.237 ±0.079 




2 


409.541 ± 0.033 


0.133 ±0.018 


2 


439.290 ± 0.016 


0.182 ±0.016 


2 


468.457 ± 0.087 




2 


498.239 ± 0.089 




2 


528.407 ± 0.374 




3 


420.000 ±0.109 
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To remedy this. lDeheuvels fc Michell (|201ll ) proposed a method specifically designed to handle these 
stars. They showed that the combined knowledge of the mean large separation {Av) of the star 
and the frequency of an avoided crossing v cross provides extremely precise estimates of the stellar 
mass and age, for a given physics. The reason for this is that both {Av) and f cr oss are monotonic 
functions of age. As a result, there is one and only one value of the mass and age for which 
the observed values of thes e quantities are simultaneously reproduced. For each considered physics, 
Deheuvels &: Michell (|201ll ) therefore suggested to fix the stellar mass and age to the values required 
by {Av) and Across) which solves the problem caused by the short timescale of avoided crossings. 

This method was tailored for subgiants, whose spectra usually contain few avoided crossings. 
The star KIC7341231 is obviously more evolved than these objects because we had to introduce 
about 15 7 modes of degree I = 1 to account for the observed oscillation spectrum of the star (see 
Sect. 13.2ft . Consecutive avoided crossings are close to each other and it is hard to isolate them. 
The method proposed bv lDeheuvels &: Michell (|201ll ) therefore had to be adapted. 



Contrary to subgiants, the 7 modes of KIC7341231 have a high enough order to follow the 
asymptotic theory. Their periods are thus approximately equally spaced by the period spacing, 
defined as 

vr f f r A by 



AH, = - 



-dr 



-1 



(4) 



where L = \Jl(l + 1) and Abv the Brunt- Vaisala frequency. Since ./Vbv keeps increasing with age 
due to the growing central density, the perio d spacing AIT/ mono t qnical ly decreases as the star 
evolves. Therefore, in the method proposed bv lDeheuvels Michell (|201ll ). the quantity f cross can 
legitimately be replaced by Alii. The observed values of {Av) and (Alii) can be used to obtain 
good first estimates of the stellar mass and age, for any given physics. This point is illustrated 
in Fig. [3 The mass and age can then be fine-tuned to reproduce at best the observables of the 
star. We note that since {Av) and (Alii) are used only to derive a first estimate of the stellar 
parameters, we can safely use the values derived for these quantities in Sect. 



4.2. Modeling of KIC7341231 

We applied the method described above to model the star KIC7341231. 



4-2.1. Properties of the models 



All the models were computed with the evolution code C ESAM2K ([M orel 1997]). We used the 

The nuclear 



OPAL 2005 equation of state and opacity tables as described i nlLebreton et al 



(2008). 



reaction rates were computed using the NACRE compilation ( Angulo et al. 19991 ). The atmosphere 
was described by Edding ton's grey law. We assumed the classical solar mixture of heavy elements of 



Grevesse fc Noelsl (|1993l ). Convection was treated using the Canuto-Goldman-Mazzitelly (CGM) 
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formalism (jCanuto et al.lll996l ). This description involves a free parameter, the mixing length, 
which is taken as a fraction cvcgm of the pressur e scale height H p . In this work, we assumed a 
value of «cgm calibrated on the Sun (a© = 0.64, ISamadi et al.ll2006l ). The effects of microscopic 
diffusion were neglected in this study. 



We used the oscillation code LOSC (jScuflaire et al.l l2008l ) to compute the mode frequencies 
of the models. It is well known that absolute mode frequencies are affected by our improper 
modeling of surface convection (see e.g. IChristensen-Dalsga ard &: Thompson! 1 199 71 ). We used the 
correction of surface effects proposed by Kjeldsen et al. ( 20081 ). which consists of adding to the 
mode frequencies a power law whose exponent is calibrated on the Sun. With our evolution code 
and treatment of convection, we found an exponent of 4.25, which was used for all the models in this 
work. We also note that some modes in the spectrum of the star are mixed and are therefore less 
sensitive to surface effects. To take this into account, the surface correction of non-radial modes 
was multiplied by a factor Q~), where Q n i corr esponds to the rat io of the mode inertia to the 



inertia of the closest radial mode, as prescribed bv lAerts et al.1 (120101) CCha p . 7). We note that the 
empirical way of correcting for surface effects proposed by Kjeldsen et al. (j2008l ) is probably not 
optimal, and efforts are currently made to take them into account in a more satisfactory way (e.g. 
Gruberbauer et al.ll2012l ). 



4.2.2. Results 

We applied the method described in Sect. 14.11 to find stellar models reproducing both the seis- 
mic and atmospheric constraints of KIC7341231. Since the metallicity of the star is rather uncer- 
tain, we tried several values of (Z/X) ranging from —1.75 to —0.75 dex. The ini tial helium content 



was e stimated from galactic evolution, by supposing a ratio A Y/AZ = 3 ± 2 (IPaee 



&: Portinari 



19981 ) and an primordial helium abundance of Y p = 0.2477 (jPeimbert et al.1 120071 ) . For each 
value, we determined the stellar mass and age for which we both have {Au) = 28.9 ± 0.2 ^uHz 
and (AIIx) = 107.1 ± 2.3 s, as was obtained for the star in Sect. 13.21 We then fine-tuned these 
parameters to match the atmospheric constraints (T e g and log g) and the frequencies of the ob- 
served modes as closely as possible. We have here used our photometric estimate of the effective 
temperature (T e g = 5470 ± 150 K, see Sect. I2.2p . but we verified that our conclusions remain 
unchanged if we consider our spectroscopic estimate of T e g instead. The agreement between the 
models and the observations was estimated using the reduced \ 2 function defined as 

1 N I /-/-jobs //-)mod\2 

where C£ bs , k = 1,N correspond to the N observables available for the star, their error bars, 
and 0™ od the corresponding values in the computed models. The number of seismic constraints is 
much larger than the number of atmospheric constraints. To avoid drowning the contribution of 
the latter in the total value of the \ 2 -> we computed a separate x 2 f° r t ne seismic observables (x 2 e is) 




378 380 382 384 386 

Frequency Qu,Hz) 



Fig. 4. — Section of the power spectrum of KIC7341231 in the neighborhood of the 384-^Hz radial 
mode (light grey area). The dark grey area corresponds to a rotationally-split I = 1 mode and the 
vertical arrows indicate peaks that are attributed to I = 2 modes (see text). 



Table 3: Parameters of the best-fit models as a function of the assumed metallicity. 



Model 


A 


B 


C 


D 


E 


[Z/X] (dex) 


-0.75 


-1. 


-1.25 


-1.5 


-1.75 


Mass (M ) 


0.880 


0.836 


0.804 


0.790 


0.770 


Age (Gyr) 


11.3 


12.2 


13.1 


13.4 


14.3 


T cS (K) 


5245 


5363 


5452 


5521 


5566 


Radius (Rq) 


2.67 


2.62 


2.59 


2.58 


2.55 


log g 


3.527 


3.520 


3.514 


3.510 


3.508 


r BCE C^*) 


0.33 


0.36 


0.39 


0.43 


0.46 


Xatm 


0.3 


0.2 


0.5 


1.0 


1.3 


x 2 - 


7.7 


5.5 


5.6 


6.8 


9.6 


" BCE: base of the convective envelope 
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Fig. 5. — Profiles of four I = 1 multiplets in the oscillation spectrum of KIC7341231. The profiles 
we obtained when assuming a varying splitting in our fit are overplotted in red. The fitted values 
of the rotational splittings are specified for each multiplet. 
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Fig. 6. — Profiles of two I = 2 multiplets in the oscillation spectrum of KIC7341231. The profiles 
we obtained when assuming a varying splitting in our fit are overplotted in purple. The fitted 
values of the rotational splittings are specified for each multiplet. 




Fig. 7. — Location in the HR diagram of the models that reproduce the observed value of Alii 
(red dashed line) and the observed value of Av (blue dotted line). The models represented here 
are computed with [Fe/H] = — 1 dex. The evolutionary tracks of several models, whose masses are 
specified in M & , were overplotted. 
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and for the atmospheric observables (Xatm)- 

The properties of the best-fit models that we obtained for each considered value of the stellar 
metallicity are given in Table [3J The observables are best reproduced for a metallicity around 
— 1 dex. The low value of Xatm an d Xseis f° r these models indicates a close agreement with both 
the atmospheric and the seismic constraints of the star, which is confirmed by Fig. El For lower 
abundances of heavy elements, the effective temperature of the models differs from the observations 
(this is all the more true when considering our spectroscopic estimate of T e g for the star) and more 
importantly, the agreement with the observed mode frequencies deteriorates. KIC7341231 is found 
to be a low-mass star (0.77 to 0.88 M@) at the bottom of the red giant branch. Its age (ranging 
from 11.3 to 14.3 Gyr) is consistent with that of an old halo star. The surface gravity of the models 
matches the seismic estimate of logg = 3.55 ± 0.03 that was obtained in Sect. [3] within 1-a error 
bars. We note that in all our best-fit models, there is an I = 2 avoided crossing around 380 /uHz 
(see Fig. [8]). This confirms the hypothesis made in Sect. 13.31 that two I = 2 mixed modes are 
detected in the oscillation spectrum of the star around this frequency. Finally, the spectroscopic 
upper limit that we obtained on vs'mi (< 1 ± 1 kms -1 , see Sect. 12. 3p combined with the estimate 
of the star's radius that we get from the models (about 2.6 R & ) enables us to set an upper limit 
on the surface rotation rate of ^ S urf < 88 ± 90 nHz (we recall that the inclination angle of the star 
has been found to be 85 ± 5°). 

From the models, we can compute the rotational kernels of the detected modes. Fig. [9^. shows 
the integrated kernels of three I = 1 modes trapped in different regions inside the star, computed 
using the models presented in Table [3] (only model E was excluded because its seismic \ 2 1S larger). 
The variations in the kernel profiles are quite small when switching from one model to another, 
which shows that the trapping of the I = 1 modes depends only weakly on the model we choose for 
the star. The integrated kernels of the two I = 2 modes for which the rotational splitting could be 
measured are plotted in Fig. [9b. For one of them (y = 439 /xHz), we reach the same conclusion as 
for the I = 1 modes. However, for the second mode (y = 409 /(xHz), the variations are larger. The 
reason for this is that this mode is in fact undergoing an avoided crossing with another / = 2 mode 
trapped mainly in the core. Unfortunately, this latter mode could not be detected. As a result, the 
trapping of the detected mode is uncertain and varies from one model to the other, which generates 
the differences that we observe between the rotational kernels of the models. We note that a longer 
data set will perhaps make it possible to detect the gravity-dominated I = 2 mode, which would 
solve this problem. The I = 1 modes are not prone to this kind of uncertainty because they can be 
detected even when they are trapped essentially in the core. 



4.3. Rotational splitting vs mode trapping 

Having access to a stellar model of the star, we were then able to relate the observed rotational 
splittings to the trapping of the modes, the idea being to search for the signature of differential 
rotation in the interior of the star. For this purpose, we introduced for each mode the weighted 
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radius r n j 

'W «» 



rK n j(r)dr 



which corresponds to the mean location of the rotational kernel inside the star. Small values of v n i 
indicate that the mode is trapped in the core and therefore more g-like, whereas if v n i is close to 
1, the mode is trapped in the envelope and close to an acoustic mode. 

In Fig. [10l we plotted the measured splittings as a function of the weighted radii r n> i, that were 
computed for each detected mode using model B as a reference model. We deliberately omitted 
the I = 2 mode around v = 409 /^Hz in this figure because it is the only mode for which the value of 
the ratio q depends significantly on the choice of the reference model, for the reason we mentioned 
in Sect. 14.2.21 We observe a very clear correlation between these quantities: the modes which have 
a g-mode behavior tend to have a larger splitting than the ones which behave as p modes. This 
clearly suggests that the core of KIC7341231 rotates faster than the envelope. 



5. Inversion of the rotational profile O(r) 

We established that the core of KIC7341231 is rotating faster than the surface. In order to 
quantify this differential rotation, we then tried to estimate the rotation profile of the star by 
inverting Eq. [21 using the observed splittings and the rotational kern els of the modes computed in 



Sect. HI We start by recalling general background on inversions (see IChristensen-Dalsgaard et al 



1990 for more details). 



We performed linear inversions, which means that for each radius tq inside the star, the inferred 
rotation profile Cl{tq) can be expressed as a combination of the rotational splittings 



M 



^( r o) = c k(r )Suj k (7) 



fc=l 

where we have used a subscript k = 1,M for the detected modes instead of their radial order n 
and degree I, for convenience. The coefficients Cfc(ro) are either determined during the process of 
the inversion or they can be calculated separately. The solution is characterized by the averaging 
kernels defined as 

M 

JC(r;r ) = Y,Ck(r )K k (r) (8) 
fc=l 

where K k (r) is the rotational kernel of the k th mode. By combining Eq. [2] and we obtain that 
^( r o) = Jo ^-( r ; r o)^( r ) dr. We thus aim at getting the averaging kernel JC(r;ro) as localized as 
possible around ro to ensure that fi(ro) is a good approximation of the true rotation rate O(ro). 
In all the cases below, the averaging kernels are required to have unit integral, so that the solution 
fj(ro) constitutes a proper average of the true rotation rate and its standard deviation is then given 
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Fig. 8. — Echelle diagram of the mode frequencies of KIC7341231 obtained in Sect. I3.3I (filled black 
circles). We overplotted the mode frequencies of the best-fit model obtained with [Z/X] = — 1 
dex (blue squares: I = modes, red triangles: I = 1 modes, purple diamonds: I = 2 modes, grey 
plusses: i = 3 modes). 
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Fig. 9. — Left: Cumulative normalized integral of the rotational kernels of three I = 1 modes that 
were detected in the spectrum of KIC7341231, computed using models A (full lines), B (dotted 
lines), C (dashed lines), and D (dash-dot lines). The black curves correspond to a mode mainly 
trapped in the envelope, the blue curves to a mode mainly trapped in the core and the red curves 
to a mode in between. Right: Same as left panel for the two I = 2 modes for which the rotational 
splitting could be determined (black: v = 409 //Hz, blue: v = 439 /iHz). 
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by 



where the ag u k are the standard deviations of the observed modes. 



M 
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\ fe=l 



(9) 



We have tried two of the most commonly used inversion techniques to invert the rotation profile 
of the star, the Regularized Least Squares (RLS) method and the Optimally Localized Averages 
(OLA) method. These methods were both successfully appli ed to estimate the internal rotation 



profile of the Sun (e.g. ISchou et al. 



1998 



Chaplin et al.lll999i ). In what follows, we have used the 



mode rotational kernels computed from model B, which is the one that offers the closest fit to the 
observations, and we have excluded from our list of measured splittings the / = 2 mode around 
409 /iHz because its rotational kernel is too model-dependent. The results that are described are 
very similar if we use models A, C or D instead of model B. 



5.1. Least Squares methods 

One approach to inverting the rotation profile of the star is to try to reproduce as closely as 
possible the observed rotational splittings by performing a Least-Squares fit to the observations. 
The objective is thus to minimize the x 2 function 

M \Su k -f R ti(r)K k (r)dr} 2 

* = E 1 -2 ~ ( 10 ) 

k=l 8u),k 

In practice, this can be done by discretizing the rotation profile £l(r) on an Appoint grid. The 
stellar radius is split into N regions delimited by the radii = ro<...<rjv = l and for j = 1,N 
we set 

Cl(r) = for Tj—i < r ^ rj (11) 

In the following, we first tried small values of N, resulting in simple profiles. Then, for larger values 
of N, the solution needs to be regularized to ensure numerical stability. 



5.1.1. Solid-body rotation profiles (N = 1 ) 

The case N = 1 corresponds to a uniform rotation rate O(r) = Sli throughout the star. By 
minimizing the \ 2 function given by Eq. [TPl we obtained an optimal value of Vl\ = 328 ± 7 nHz. 
The agreement with the observed data is very poor (blue circles in Fig. [10]) . yielding a reduced x 2 
of about 17. This clearly shows that we can reject the hypothesis of a solid-body rotation inside 
the star. Fig. [10] shows that if we assume a uniform rotation, the modes trapped in the core 
are expected to have slightly smaller splittings than the modes trapped in the envelope — this is 
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caused by the fact that the kernels of p modes have approximately unit integral, while the integral 
of g-mode kernels is smaller than one. In the observations, it is the contrary, which confirms that 
the core rotates faster than the envelope. 

5.1.2. Two- zone models (N = 2) 

For N = 2, the star is split into two uniformly rotating regions separated at a radius with 
fj(r) = f^i toward the core and Cl(r) = O2 toward the surface. 

We first separated these zones at the interface tce between the convective envelope and the 
radiative core. We note that based on the model we consider, this value of tce varies (see Table 
[3]). However, the results that we obtained for the rotation rates vary little. By minimizing the x 2 
function, we obtained {l± = 696 ± 24 nHz and SI2 = 51 ± 19 nHz. The agreement between the 
theoretical splittings and the observed ones is much better than with the solid-body profile, as can 
be seen in Fig. \W\ The reduced x 2 ls now about 1.6, which indicates a fairly good agreement 
with the observations. To estimate how robust this result is, we repeated the same procedure 18 
times, rejecting points one by one in our list of measured rotational splittings. Fig. [TT] shows 
that the values of £li and 0,2 remain quite stable in all cases. We then calculated the averaging 
kernels corresponding to the solution that we obtained. Instead of plotting the averaging kernels 
themselves, we chose to represent in Fig. [12] the cumulative integral of their modulus in order to 
estimate more efficiently to which region of the star they are most sensitive (we recall that unlike 
the rotational kernels of the modes, the averaging kernels are not necessarily positive functions). 
Apart from a small contribution from the surface, the core kernel is quite well localized. It is 
especially sensitive to the innermost 2% of the star in radius. On the contrary, the envelope kernel 
appears to be severely contaminated by the core. 

We also tried to consider the intermediate radius r\ as a free parameter instead of fixing it 
to rcE- In this case, the minimization is not as trivial since Vx 2 is not a linear function of the 
rotation rates at the grid points. However, the simplicity of the profile allowed us to vary r\ from 
to 1 and each time determine the optimum values of f2i and fl2- Fig. [13] shows the variations of 
the reduced x 2 that we obtained as a function of r\ . The first remark is that the minimum x 2 1S as 
low as 0.93, significantly smaller than the one obtained with n = tce- The corresponding rotation 
profile is very strange. We obtained n = 0.985, VLi = 730 ± 27 nHz and 2 = -2513 ± 178 nHz, 
i.e. a flat rotation in the largest part of the star and a thin layer at the surface spinning fast in the 
opposite direction. 

However, there are several indications showing that this solution should be discarded. First, 
the x 2 function has an oscillatory behavior and thus many secondary minima throughout the star's 
interior. Secondly, knowing that the radius of the star is about 2.6 Rq, a surface rotation rate of 
1^2 = —2513 nHz would require a surface velocity of about 29 kms -1 , clearly inconsistent with the 
spectroscopic upper limit of v sini < 1 ± 1 kms" 1 . Finally, we performed a simulation to show that 
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Fig. 10. — Rotational splittings of the modes as a function of the weighted radii r n i given by Eq. 
[6J The empty symbols (1 = 1: circles, 1 = 2: squares) correspond to the observations with 1-a error 
bars, the filled blue symbols to an optimal solid-body rotation model and the filled red symbols to 
an optimal two-zone model. 





660 670 680 690 700 710 720 730 
Ci, (nHz) 



Fig. 11. — Values of Oi and obtained when rejecting points one by one in the rotational splittings 
list (black crosses). The red cross indicates the result obtained when considering all the rotational 
splittings. 
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Fig. 12. — Integrated modulus of the averaging kernels of the core (dashed lines) and the surface 
(full lines) obtained by describing the rotation profile as a two-zone model (black), by using the 
RLS inversion method with a 100-point grid (blue) or the OLA inversion method with the same 
grid (red). For more clarity, the integrals were normalized to their surface value. 
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Fig. 13. — Left: Variations in the reduced x 2 value as a function of the radius T\. The red circle 
indicates the minimum of the function. Right: Optimal values of the core rotation Vt\ (red) and 
the surface rotation Q2 (blue) as a function of r%. The dashed lines indicate the error bars. 
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the minimum at r% = 0.985 is spurious. We computed theoretical rotational splittings using the 
optimal rotation profile that had been obtained when setting n = tce (=0.36 for model B). We 
added to these theoretical values a Gaussian noise whose width corresponds to the error bars of the 
observed splittings. We then tried to recover the input value of r\ by applying the same procedure 
as before. We performed several iterations of this procedure, which all led to a x 2 function very 
similar to the one that we obtained from the data. They indeed have an oscillatory behavior and 
their minimum can be found anywhere between r = and r = 1 instead of recovering the input 
value rcE = 0.36. Fig. [TH shows the example of an iteration for which the minimum of the x 2 
function is close to the surface. All this shows that we can ignore this strange counter-rotating 
profile. 



5.1.3. Regularized Least Squares method 

We then explored the case of larger values of N. In this case, the finite number of measured 
splittings for the star is obviously not sufficient to reconstruct the whole function fi(r), which 
makes any solution non unique. Besides, the solutions are very sensitive to both the discretization 
and the measurement errors (Tguk- As a result, the problem is ill-conditioned and needs to be 
regularized. The Regularized Least Squares method consists in adding to the classical x 2 function 
to be minimized a well chosen regularization function. 



To apply the RLS method to our star, we followed the procedure prescribed bv lChristensen-Dalsgaard et al 



()1990l ). The regularization term was chosen as a smoothness constraint on the solution. We thus 
minimized the function 

J = X 2 + l*uF(il) (12) 

where F(£l) is the regularization function, which was taken as the norm of the second derivative of 
fi(r), i.e. F(fl) = \\£l"(r)\\ 2 . The factor \x is a trade-off parameter between the smoothness of the 
solution and the minimization of the x 2 function. For more convenience, we added a normalization 

factor , defined u ^ (t11\ Ju\\ The minimization of the function J can be written as 



a linear problem (jChristensen-Dalsgaard et al.lll990l ). which can readily be solved to estimate the 



optimal rotation profile. 

To determine a satisfactory value of the trade-off parameter fj,, we generated artificial rotation 
profiles fn(r), from which we computed theoretical rotational splittings. We then applied the method 
described above with different values of the regularization parameter [i to try to recover the input 
rotation profile. We found that a value of /i « 100 provides a good compromise. Smaller values of 
H introduce spurious large-amplitude oscillations in the recovered profile, while larger values of ji 
lead to an over-regularization and the obtained profile corresponds to the one that minimizes F(Q) 
(a straight line if the norm of the second derivative of O(r) is taken as a smoothness condition). 
We therefore used n = 100 when inverting Cl(r) from the observed splittings. 

Fig. [15] shows the rotation profile that we obtained by discretizing f2(r) over a 100-point grid. 
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Fig. 14. — Same as the left panel of Fig. [131 but for a simulation (see text). 
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Fig. 15. — Left: Inverted rotation profile obtained with the RLS method (100 points, regularization 
with the second derivative and \x = 100). The dashed lines indicate 1-a error-bars. We note that 
apart from in the core, the averaging kernels are ill-localized and the inverted rotation rates are 
not reliable (see text). Right: Same as left panel using the OLA method. 
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We found a rotation rate in the core of 695 ± 29 nHz, very similar to the one derived with two- 
zone rotation models. A rotation rate of —272 ± 137 nHz was obtained at the surface of the star. 
The agreement between the rotational splittings obtained with the inverted rotation profile and 
the observed ones is only slightly better than the one we reached with a two-zone model in Sect. 
15.1.21 Useful information about the quality of the inverted profile can be obtained by computing 
the corresponding a veraging kernels. Details on how t hese kernels are calculated with the RLS 
method are given by IChristensen-Dalsgaard et al.l ()1990l ) . An inspection of these averaging kernels 
shows that only the one located in the core has a satisfactory shape. It is in fact very similar to 
the core kernel obtained with the two-zone model. It is well localized in the innermost 1.4% of the 
stellar radius (17% in mass) and has a small contribution from the surface (see Fig. [T6|) . On the 
contrary, all the other kernels have large leakage from both the core and the surface. The integrated 
modulus of the surface kernel, for instance, has a large contribution from the core (see Fig. [T2]) . 
Only the estimate of the rotation rate obtained for the core is thus reliable. In particular, the slight 
counter-rotation that we obtained at the surface cannot be regarded as significant. 



5.2. Optimally Localized Averages method 

Instead of seeking to reproduce at best the observed rotational splittings like the RLS method, 
the Optimally Localized Averages (OLA) method consists in calculating localized averages of the 
true rotation profile in different regions of the star. For this purpose, the method builds for each 
point rj inside the star a linear combination of the mode kernels ^2k c k( r j)^-k( r ) such that the 
resulting averaging kernel /C(r;r 3 ) is as close as possible to a Dirac function centered on rj. 



Following iBackus Gilbertl (119681 ). we searched for the coefficients Ck{rj) by minimizing the 
function 

[-R M 
J = 12/ IC^rtfir-rrfdr + fiJ^ickir^as^} 2 (13) 

with the constraint that J JC(r;rj) dr = 1, for each point rj; /x is a trade-off parameter between 
resolution of the averaging kernels and error magnification. 

In our case, it is very hard to obtain localized averaging kernels, so we chose to set \i = 0. 
Even this way, most averaging kernels suffer from very large leakage from other parts of the star 
(especially from the core and the surface). The core kernel is an exception (see Fig. [TBI) . Fig. [12] 
shows that the leakage from the surface in the core kernel that was observed with the RLS method 
has now been cancelled. If we assume that the rotation profile of the star varies smoothly in the 
core, then the core kernel can be well approximated by a step function H(r) between and a small 
radius r ~ 0.014 R (see Fig. [T6|) . This means that 

^,f"a„ (r)Jr= w» (14) 

Jo -ft r o 
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The quantity 0(0) is therefore a very good approximation of the average of O(r) in the innermost 
1.4% in radius, which corresponds to 17% in mass, of the star. We obtain 710 ± 51 nHz for this 
average, which is consistent with the result of the RLS method. 

Fig. fl2l shows the integrated value of the surface averaging kernel |/C(r;l)| along the stellar 
radius. The contribution from the core to this kernel has been significantly decreased compared to 
the RLS method. However, there remains a small leakage that can have significant consequences 
because we know that the core rotation is larger than the envelope rotation. The estimate of the 
surface rotation (150 ± 19 nHz) is therefore certainly overestimated and only gives an upper limit 
of the surface rotation. 

6. Conclusion 

In this paper, we obtained a precise seismic determination of the rotation rate in the core of the 
early red giant KIC7341231 and we proved that it spins at least five times faster than the surface. 

KIC7341231 is a low-mass evolved star, which is located at the bottom of the red giant branch. 
Solar-like oscillations have been detected in the oscillation spectrum of the star, derived from one 
year of Kepler observations (quarters Q5-6-7-8). Due to the evolution stage of the star, many of 
its non-radial modes have a mixed nature, which means that they behave both as g modes in the 
core and as p modes in the envelope. We found that many of these mixed modes are very clearly 
split by stellar rotation and we therefore set out to probe the rotation profile of the star. 

We performed a seismic analysis of the oscillation spectrum of the star and were able to 
determine precisely the rotational splittings of 19 I = 1 and I = 2 modes. They were found to 
range from 0.13 to 0.41 /iHz, with error bars of 0.03 fjMz on average, thereby suggesting that the 
interior of the star is differentially rotating in radius. We then found a stellar model reproducing 
very well both the atmospheric and the seismic properties of the star. We used this model to study 
the relation between the observed rotational splittings and the regions in the star where the modes 
are trapped. We found a clear correlation between these quantities that unambiguously indicated 
that the core rotates faster than the envelope in KIC7341231. 

Finally, we performed inversions of the rotation profile of the star, using the observed splittings 
and the rotational kernels of our best-fit models. We used both the RLS (Regularized Least Square) 
and the OLA (Optimally Localized Averages) methods and obtained the following results: 

• We were able to determine a very robust and precise estimate of the core rotation of the 
star. All the methods that we used (RLS, OLA) provided a core rotation rate consistent with 
f2 c = 710 ± 51 nHz within 1-a error bars. Besides, we obtained similar values when using 
the rotational kernels of other models of the star computed with different metallicities, so 
this result seems to be only weakly model-dependent. We showed, using the OLA method, 
that this rotation rate in fact corresponds to a very good approximation of the average of 
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f2(r) in the innermost 1.4% of the stellar radius. It is ironic that the core rotation rate of 
KIC7341231 could be measured while the solar core rotation is still uncertain for r < 0.2Rq 



(jChaplin et alJll99a ) 



We obtained an upper limit for the surface rotation of f2 s < 150 ± 19 nHz. This enabled us 
to establish that the core rotates at least five times faster than the surface in this star. 



We note that in this study we essentially focused on the rotational splittings of I = 1 modes. 
This is partly caused by the SNR of modes of degree I ^ 2, which in most cases remains too low to 
reliably determine their rotational splittings. This problem should be at least partially solved by 
the growing data set from the Kepler spacecraft (this star will continue to be on the short-cadence 
target list at least through Q12 and hopefully for the remainder of the mission). However, we also 
showed that the profiles of certain / = 2 mixed modes split by rotation significantly differ from the 
expected one and it is our opinion that some theoretical work still remains to be done to better 
understand the effects of rotation on mixed modes of degree / ^ 2. 

We are entering a new era in the study of the transport of angular momentum in stars because 
we now have access to observational constraints on the internal rotation profiles of stars, which were 
longed for since a very long time. We will very likely find among the Kepler targets other subgiants 
and red giants whose internal rotation can be inferred by interpreting the rotational splittings of 
mixed modes. For instance, if we assume that the rotation is nearly rigid at the end of the main 
sequence (as it is in the Sun), then the differential rotation observed in subgiants and red giants is 
entirely caused by the core contraction in the post main sequence stage. The ratio between the core 
rotation and the surface rotation at different luminosities along the giant branch should bring very 
valuable constraints on the timescale of the exchange of angular momentum between the core and 
the envelope. This should help us determine which mechanisms of angular momentum transport 
dominate and how efficient they are. 
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A. Modeling of the background in the Power Spectral Density 



The background B(v) is modeled as the sum of a white noise B w describing the photon noise 
(which dominates the spectrum at high frequency ) and a Harvey profile Bg{v) modeling the gran- 
ulation spectrum. As prescribed by lHarvevI (|19851 ). we assumed for B g {v) a function of the type 

4(7^ TV 

BJv) = , g g r (Al) 

gV ' l + (2^T g ) a s V ^ 



This preprint was prepared with the AAS IATjjX macros v5.2. 
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where r g and a s cor respond to the characteristic timescale and amplitude of the granulation (see 



Mathur et al.ll2011bl for a discussion on methodologies to fit the background). The exponent a g 
corresponds to the slope of the power law at high frequency. In most analyses of solar-like oscillation 
spectra, a function B{y) is fitted to the background prior to the fit of the p-mode component. This 
has the advantage of stabilizing the latter fit and it has been shown several times that this has neg- 
ligible influence on the results. While fitting the background component, the p-mode contribution 
is modeled as a Gaussian function. We used an MLE (Maximum Likelihood Estimation) method 
to fit the background, taking into account the fact that the bins of the PSD follow a 2-degree-of- 
freedom \ 2 distribution. The fitted background is shown in Fig. [5] and the obtained parameters 
are given in Table HI As a by-product, we obtain an estimate of the frequency of the maximum 

signal l^max 

= 406 db 3 //Hz, which corresponds to the centroid of the Gaussian used to describe the 
stellar pulsation contribution in the power spectrum. 



Table 4: Fitted parameters of the background model B(v) (see text). 



cig (10 3 ppm) r g (s) a g B w (ppm 2 //iHz) 

99 ±5 1096 ± 59 2.17 ±0.05 3.32 ± 0.04 
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Fig. 16. — Left: Averaging kernel in the center of the star obtained with the RLS method (black 
curve) and the OLA method (red curve). Right: Zoom of the left panel in the core. The blue 
dashed line corresponds to the step function H(r) (see text). 



